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Abstract 

We present an analysis of the composition of inclusive semileptonic B meson decays using 9.4 fb~^ 
of e+e" data taken with the CLEO detector at the T(4iS') resonance. In addition to measuring the 
charged lepton kinematics, the neutrino 4-vector is inferred using the hermeticity of the detector. 
We perform a maximum likelihood fit over the full three-dimensional differential decay distribution 
for the fractional contributions from the B XJiy processes with Xc = D, D*, D**, and non- 
resonant Xc, and the process B X^lv. Prom the fit results we extract = (4.05 zt 0.18 =b 
0.58 lb 0.25 lb 0.21 lb 0.56) x 10~^ where the errors are statistical, detector systematics , B — > X^lv 
model dependence, B — > Xylu model dependence, and theoretical uncertainty respectively. 



*Submitted to the 31®* International Conference on High Energy Physics, July 2002, Amsterdam 



2 



I. INTRODUCTION 



The CKM matrix element \Vub\ is the couphng between the bottom and up quarks, and 
is therefore a fundamental parameter in the Standard Model. Its value is relevant to studies 
of flavor-changing currents, including the rates of B and mixing and CP violation in the 
mixing and decays of B hadrons. 

Recent theoretical progress in charmless semileptonic B decays using Heavy Quark Ef- 
fective Theory (HQET) and the Operator Product Expansion (OPE) ^ ^, ^ has made 
possible new levels of precision in the measurement of |V^fe|[|l- The calculations relate \Vub\ 
to the partial branching fraction of B Xulv in a restricted region of phase space in which 
B — >■ Xclv does not contribute. The region of phase space most accessible to experiment 
is the lepton energy endpoint, Ei > {m\ — m'j^)/2mB- This region was used in the first 
measurement of B ^ XJ,v establishing a non-zero value of Unfortunately the 

OPE calculation is not valid in the very restricted region of phase space delineated by this 
lepton energy cut. A recent CLEO measurement has used a spectral function measured in 
B — > to obtain the normalization factor which relates the measured partial branching 
fraction in the lepton endpoint region to |Vu6|||^. Another method is to replace the lepton 
energy cut with a cut on the invariant mass of the lepton pair0, |^, g^, or on the mass of 
the recoiling hadronic system, Mx, to exclude B — > Xclv events. Recent theoretical work 
suggests that a combination of cuts on both and Mx yields the smallest uncertainty^. 

In this paper we report an analysis of i? — > Xlv decays in which both the charged lepton 
and the neutrino are reconstructed. The neutrino is reconstructed using the approximate 
hermeticity of the CLEO II and CLEO II. V detectors and the well known initial state of 
the e^e~ system produced by the Cornell Electron Storage Ring (CESR). The kinematic 
variable is then calculated directly, and the Mx kinematic variable can be calculated if the 
B momentum, ps, is known. 

Ml = Ml + q^- 2EbeaUEe + E^) + 2\pB\\Piv\ cosOB-eu, 

where pei, is the sum of the charged lepton and neutrino momenta and 6b £u is the angle 
between pi^, and the B momentum direction. Since the B mesons are the daughters of an 
T(4S') produced at rest, the magnitude of the B momentum is known and small, however its 
direction is unmeasured. The last term in the M^ equation depends on the B momentum 
direction, and is small, unmeasured, and neglected in this analysis. From the lepton pair 
sample, the differential decay rate is measured as a function of the reconstructed quantities 
g^, Mx, and Ef. The experimental resolution on the neutrino four-momentum is poor, with 
a narrow core of approximately 120 MeV and broad tail of over-estimation of the neutrino 
energy which extends up to 1.5 GeV. The resulting resolution on and Mx is also poor, so 
it is not possible to cleanly isolate B X„/z/ events with cuts on these variables. Instead, 
we fit this observed distribution to a sum of models for the different hadronic final states: 
D, D*, D**, Xc non-resonant, and Xu- From the fit results we infer the partial branching 
fraction over a region of phase space with high and low Mx- These partial branching 
fractions are related to \ Vub\ by theoretical calculations whose uncertainty has been estimated 

i- 

The data were accumulated with two configurations of the CLEO detector |^, CLEO II 
and CLEO II. V, with an integrated luminosity of 9.4 fb~^ taken on the T(45') resonance 
and an additional 4.5 fb~^ taken off resonance, just below the BB production threshold. 
Both configurations cover 95% of the Air solid angle with drift chambers and a cesium iodide 
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calorimeter. In addition there are muon chambers with measurements made at 3, 5, and 7 
hadronic interaction lengths of iron and a time of flight system which augments the particle 
ID information from specific ionization [dE/dx). In the the CLEO II configuration, there 
were three concentric drift chambers filled with a mixture of argon and ethane. In the 
CLEO II. V detector, the innermost tracking chamber was replaced with a three layer silicon 
detector and the drift chamber gas was changed to a mixture of helium and propane. 

II. NEUTRINO RECONSTRUCTION AND EVENT SELECTION 

Events are selected to have an electron or muon with momentum greater than 1 GeV/c 
and a well reconstructed neutrino. Additional cuts are used to suppress background events 
from the e^e" — > qq continuum. 

The leptons are selected to fall within the barrel region of the detector (|cos^^| < .71, 
where 9 is the angle between the lepton momentum and the beam axis). Electrons are 
selected based on E/p, dE/dx, and time of flight information combined using a likelihood 
technique. Muons are identified by requiring a penetration of 3 interaction lengths of iron 
for muons less than 1.5 GeV/c and 5 interaction lengths for muons greater than 1.5 GeV/c. 
The absolute lepton identification efficiencies arc calculated by embedding leptons from 
radiative QED events into hadronic events. The rate at which pions and kaons fake leptons 
is measured by reconstructing Kg — > tt+tt", D'^ ir~7r+, and — > X+tt" without using 
particle identification and then checking daughter particle lepton identification signatures. 

Neutrinos are reconstructed by subtracting all observed track and neutral shower four- 
momenta from the four-momentum of the e'^e~ initial state which is at rest in the laboratory, 

Pv Pe+e~ P visible' 

The errors made in this assumption are due to particles lost through inefficiency or accep- 
tance, fake tracks and showers, and other undetected particles such as second neutrinos, 
X^'s, or neutrons. Several cuts are made to select events in which these effects are reduced 
and the neutrino four-momentum resolution is correspondingly enhanced. Since extra neu- 
trinos are correlated with extra leptons, events with a lepton beyond the signal lepton in 
the event are rejected. 

The primary source of fake tracks is from charged particles which do not have sufficient 
transverse momentum, pt, to reach the calorimeter and therefore curl multiple times in 
the tracking chambers. The portions of the track after the initial outbound section may 
accidentally be called a second track. Criteria have been developed to identify such errors 
and make a best estimate of the actual charged particles in the event. Events for which the 
total charge of the tracks is not zero are removed, reducing the effect of lost or fake tracks. 
Showers in the calorimeter associated with tracks in the drift chamber are removed, so as not 
to double count their energy. The association between tracks and showers has been refined 
to take into account secondary hadronic showers which can appear to be showers isolated 
from tracks. Finally a requirement is made that the ratio of the reconstructed neutrino 
four-momentum divided by twice the reconstructed neutrino energy be less than 0.35 GeV, 

"^""'^ < .35 GeV. This quantity is constructed to be proportional to the energy of a lost 
or fake particle. The reconstructed neutrino's energy is assigned to be the magnitude of the 
momentum, because the momentum is not dependent on the particle identification of the 
tracks and so has a better resolution than the direct energy measurement. 
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Continuum events are suppressed by a combination of event shape and orientation vari- 
ables which exploit the fact that continuum events tend to be jet-like and aligned with the 
beam axis, whereas BB events are more spherical and their orientation is uniformly dis- 
tributed in the detector. The second Fox- Wolfram moment R2, of the energy flow in the 
event is required to be less than 0.4. In addition a neural network is used to combine R2, the 
angle between the lepton and the thrust axis, the angle between the lepton and the beam 
axis, and the fraction of the energy lying in 9 separate cones around the lepton direction, 
which cover the full the 4n solid angle. The R2 cut is more than 99% and 95% efficient for 
B Xclv and B —>■ Xulv respectively, while removing 60% of the continuum events. The 
neural net cut removes an additional 73% of the continuum background while keeping 92% 
and 94% of the B — >■ Xclv and B —>■ Xulv respectively. 

After all cuts we observe 42333 events from CLEO II and 81368 events from CLEO II. V. 
The overall efficiency varies from 1.5% for B XJu non- resonant to 4.2% for B —>■ XJ,v. 



III. FITTING FOR COMPOSITION 

The full three dimensional differential decay rate distribution as a function of the recon- 
structed quantities g^, Mx, and Ei is fit for the contributions of the signals and backgrounds. 
The signals are B —>■ Dlu, B D*lv, B — > D**lv, B —>■ XJu non-resonant, and B — > XJ,i^. 
The backgrounds are events where the B decay candidate is faked by a lepton from the pro- 
cess h c Xiu, a lepton from the e~^e~ —>■ qq continuum, or a fake lepton. We have the 
freedom to choose the kinematic variables used in the fit. We use q^ /{E^ + Ei,Y, Mx, and 
cos 9w£- The helicity angle of the virtual W, cos 9w£, is defined to be the angle between the 
lepton momentum in the W frame and the W in the lab frame. The fit variables have been 
chosen to minimize complexity of the phase space boundaries and cover the same kinematic 
space as g^, Ei, and E,^. 

We perform a binned maximum-likelihood fit where the probability distribution functions 
(PDFs) are constructed from weighted Monte Carlo or data events. The fit uses electrons 
and muons simultaneously, with a separate set of PDFs for the electrons and muons. The 
likelihood is implemented to take into account the PDF statistics using the method described 
in reference 0. The projections of the reconstructed quantities g^, Mx, and cos 9we of the 
various B — > Xlu modes are shown in Figure ^ Projections of the fit result are shown in 
Figure ^. Projections restricted to the region of enhanced B — Xuliy sensitivity are shown in 
Figure ^. It should be noted that the correlations between the the three variables shown in 
the projections are fully included in the fit, and provide considerable power in distinguishing 
the contributions of the various components. 

The B —>■ Xlv modes, secondary leptons and real leptons from continuum are simulated 
with CLEO's GEANT based Monte Carlo. The B Dlu and B D**lu modes are 
modeled with ISGW2 |lOl. The B ^ D*lu mode is modeled with HQET using CLEO's 



previous measurement of the form factors The Xc non-resonant modes are modeled 



with a combination of the Goity and Roberts model|jT2| and a crude model of B ^ AcXli/. 
For an assessment of the B XJi' model dependence, the B — > Dlu and B D*li> form 
factors are varied within their experimentally allowed ranges, and the B D**lv masses are 
varied in an ad-hoc manner. Because very little is known about the X^ non-resonant mass 
spectrum, it is reweighted drastically in the study of the B — >• Xclv model dependence. In a 
further study, each X^. final state is constrained to a value away from the nominal fit result. 
The results of the variations are summed in quadrature. The modeling of the B XJ,u 
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FIG. 1: Distributions of the various B Xlv modes as functions of the reconstructed quantities 
(left to right) q^, Mj^, and cosOwi- The modes are B — > Xulv (sohd), B Dlu (short dash), 
B D*lv (dots), B D**lu (dot-dash), and B XJiy non-resonant (long dash). 




FIG. 2: Projections of the fit results in in the reconstructed quantities (left to right) q^, M^, and 
cos 0\Yi- The solid histogram on the bottom is the sum of the backgrounds: continuum, secondary 
leptons, fake leptons. The hatched histogram is the sum of the B XJi' modes: B Dli/, 
B — > D*lv, B —J- D**lv, and B XJiy non-resonant. The top histogram, barely visible, is the 
B — > Xuli' component. 



mode and the resulting model dependence is discussed in a later section. 

The fake leptons are modeled by taking a sample of tracks from data that are classified as 
71, K, or (electrons and protons are not a significant source of fake leptons) and unfolding 
their spectra to get the true spectra of vr, K and /x, which are then multiplied by the 
measured fake rates. This models fake leptons from both BB and e+e~ qq continuum. 
The real leptons from continuum are modeled with continuum Monte Carlo which has been 
tuned to replicate the appropriate charm spectra. Charm is the source of most leptons from 
continuum. The models of both continuum and fakes have been validated and constrained 
by a comparison with the 4.5 fb~^ of off-resonance data. The secondary leptons are modeled 
with a convolution of CLEO's measured charm spectra from BB events and DELCO results 



on the inclusive semileptonic decays of charm|13 



6 




FIG. 3: Projections of the fit results in reconstructed quantities g^, M'^, and cosOwe with cuts 
of q"^ > 11.0 GeV'^/c^ for the Mj. projection and Mj. < 2.25 GeV^/c^ for the . projection. 
The sohd histogram on the bottom is the sum of the backgrounds: continuum, secondary leptons, 
fake leptons. The hatched histogram is the sum of the B — > Xcli' modes: B — > Dlu, B — > D*lu, 
B D**lv, and B Xclv non-resonant. The top histogram, now clearly visible, is the B — > Xulv 
component. 



The method of neutrino reconstruction adds a large amount of kinematic information to 
each event, however it also adds significant potential for systematic errors. The resolution 
on the neutrino kinematics is affected by the models of both the signal, the other B in the 
event, and the detector response. The GEANT Monte Carlo does not perfectly simulate the 
track and shower efficiencies and fake rates, nor are B decays well enough understood that 
the inclusive particle distributions are well known. For this analysis we employ a method 
of reweighting in order to quantify the effects of these uncertainties on our results. For 
example to study effect of the tracking efficiency uncertainty on result, the Monte Carlo 
events in which tracks are lost are given a higher or lower weight in constructing the PDFs 
and the fit is repeated. The scale of the variation made can in general be constrained by 
direct measurements of the quantity being varied. One important example is that having a 
K^j^ in the event adds a tail to the neutrino resolution. The inclusive spectrum in BB 
events has been measured and can be used to constrain the spectrum. 



IV. B XJv MODEL DEPENDENCE 



We use calculations based on HQET and the OPE to extract \Vuh\ from the B Xyly 
rate in restricted regions of phase space. This is a well controlled expansion and the the- 
oretical uncertainties have been assessed by several authors 0, ||, However, the fit's 
region of sensitivity does not coincide with these regions, and, as previously mentioned, it 
is not possible to make cuts to isolate these regions, because of the poor resolution on the 
neutrino four-vector and hence and Mx- In order to calculate \ Vub\i we first make a model 
dependent inference of the partial branching fraction in a region, and then apply the HQET 
and OPE calculations. This prescription is designed to minimize the reliance on models and 
instead rely on the controlled expansion in HQET and OPE calculations. 

B — >■ Xulv models are used for two purposes in this analysis. The first is in the simulation 
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of the efficiency and resolution of the detector. This cannot be done with the HQET and 
OPE calculations, because they do not predict specific hadronic final states. The second 
use of models is to extrapolate and interpolate between the regions of high sensitivity to 
B XJ,v and regions of interest for extracting \Vuh\- The fit is primarily sensitive to low 
M'j^ and high g^, but within this region the sensitivity is strongly biased toward the lepton 
energy end-point, hence the models are relied on to extrapolate to the full range of lepton 
energies. 

The result of the fit is the fraction oi B ^ Xulv events in the sample. This is converted 
into a branching fraction using the efficiency determined by Monte Carlo simulation and 
number of BBs in the sample. The B Xulv portion of the ffi is driven by the region 
of phase space where the B XJiy contribution is not overwhelming, however the fit 
extrapolates from this region to the full space using the models. This results in a very large 
model dependence for the total rate. The effect can be reduced by calculating for each model 
a fraction, f^°^f, of events in a region of enhanced B XJi^ sensitivity. The fractions, 
/region' Calculated from the model using the true and Mj., not the reconstructed 
and M^. The product of the branching fraction from the fit result for a model and f^g^^ 
is the inferred partial branching fraction of i? ^ X^li^ for the model in region selected, 

A fomodel t2t D > V 7,,\model /"model 

-^^region = XJl^) X /,egion • 

^'^regkm significantly less dependence on the model than the inferred full branching 
fraction, B{B however it still involves a model dependent combination of 

interpolation and extrapolation from the actual region of measurement. 

The models used range across extremes from the ISGW2 model to a model with only 
non-resonant modes. The ISGW2 model contains only low mass hadronic resonances and no 
non-resonant or high mass hadronic systems. The purely non-resonant model implements the 
HQET and OPE prediction for the hadronic mass spectrum based on the B Xg'y spectral 
function, where the hadronic system is decayed via JETSET. The true spectrum is likely 
somewhere between the extremes of all high mass non-resonant and all low mass resonances. 
Two permutations of the all non-resonant model correspond to changing the spectral func- 
tion parameters, A and Ai within the bounds of the CLEO measurement Another 
model combines ISGW2 and some non-resonant events. In addition, the ISGW2 model is 
reweighted to have a higher or lower cos 6we, corresponding to more or less rate in the lepton 
end-point, or reweighted to have a harder or softer spectrum. Table | shows the results 
for the total rate and the inferred partial branching fractions for five different regions. The 
central value of AiSregion, reported in the first column of Table is the center of the range 
of models, and the model dependence uncertainty is assigned to cover the full range. The fit 
results for the branching fractions of the final states containing charm agree with previous 
CLEO measurements and will be reported in a subsequent paper. 

In the context of HQET and the OPE another fraction, f^^ion^ relates the partial branch- 
ing fraction to \Vub\ 



\Vub 



3.07 ±0.12 X 10"^ 



1/2 

Aigregion l-6pS 
.001 X /™ TB 



The fractions, f^gj^^, for the five regions used have been calculated with an evaluation of 
the theoretical uncertainty by Bauer et al. in reference [^. Some of the regions have also 
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been evaluated in reference^. The \Vub\ results that correspond to the five inferred partial 
branching fractions with their respective model dependences are shown in Table |T[ We 
choose the > 11 GeVVc^ < 1.5 GeV/c^ results central value because it has the 
least model dependence. 



TABLE I: Dependence on the B — > Xulv model of the inferred total branching fraction, B{B — > 
Xulv), and the inferred partial branching fractions, AyB^gj^'jJ, for five regions. Entries are in units 
of 10~^ 







ISGW2 & 


All Non-Resonant 






Higher 


Lower 




ISGW2 


Non-Res 


Nominal 


Low A,Ai 


High A,Ai 


Harder 


Softer ij2 


cos 9\nri 


cos 9yYi 


Total B{B ^ Xnlu) 


0.996 


1.211 


1.808 


2.289 


1.620 


0.803 


1.159 


0.956 


1.039 


g2 > 6 GeV2,Mx < Md 


0.629 


0.686 


0.835 


0.859 


0.811 


0.582 


0.655 


0.606 


0.653 


g2 > 8 GeY^,Mx < 1.7 GeV 


0.497 


0.528 


0.609 


0.617 


0.598 


0.482 


0.494 


0.480 


0.516 


g2 > 11 GeV2, Mx < 1.5 GeV 


0.308 


0.317 


0.348 


0.347 


0.346 


0.321 


0.283 


0.297 


0.319 


g2 > (Mb-Md-? 


0.283 


0.293 


0.329 


0.338 


0.318 


0.298 


0.257 


0.273 


0.293 


<?2 > (Ms - Mo? 


0.331 


0.344 


0.393 


0.410 


0.377 


0.342 


0.307 


0.319 


0.343 



TABLE IL Inferred partial branching fraction by region, A^Brcgion) with model dependence in 
the first column. \Vub\ calculated from the inferred partial branching fractions with all errors in 
the second column. The errors on are statistical, detector, B —> XJv model dependence, 
B Xulv model dependence, and theoretical uncertainty respectively. Please note, these results 
do not come from fits in the restricted regions. 



Model Error 



\Vuh\± Stat ± Detector ±6— >c±f)^-u± Theory 



(3.82 ± 0.17 ± 0.55 ± 0.23 ± 0.92 ± 0.12) x 10"^ 
(3.84 ± 0.17 ± 0.55 ± 0.23 ± 0.37 ± 0.31) x IQ-^ 
(3.98 ± 0.18 ± 0.57 ± 0.24 ± 0.25 ± 0.38) x 10"^ 
(4.05 ± 0.18 ± 0.58 ± 0.25 ± 0.21 ± 0.56) x IQ-^ 
(4.07 ± 0.18 ± 0.58 ± 0.25 ± 0.28 ± 0.62) x 10"^ 
(4.05 ill 0.18 ill 0.58 ± 0.25 ± 0.29 ± 0.54) x IQ-^ 



(1.546 ±0.743) X 10"^ 
(0.720 ±0.138) x 10-3 
(0.548 ± 0.069) X 10"^ 
(0.315 ±0.032) X 10-3 
(0.297 ±0.041) X 10-3 
(0.359 ± 0.052) X 10-3 



Total B(B Xulu) 
g2 > 6 GeV2,Mx < Mo 
i}2 > 8 GeV2,Mx < 1.7 GeV 
> 11 GeV^, Mx < 1.5 GeV 
g2 > (Mb-Md' ? 
g2 > (Mb-Md? 



V. CONCLUSION 

We have made a preliminary measurement the CKM parameter \Vub\ with high precision, 
but some model dependence, and obtain 

Vub = (4.05 ±0.18 ±0.58 ±0.25 ±0.21 ±0.56) X 10"^ 

where the errors are statistical, detector systematics, B XJu model dependence, 
B — > Xy/z/ model dependence, and theoretical uncertainty respectively. The result is con- 
sistent with the CLEO result using the lepton end-point measurement 0. There is some 
statistical correlation with the lepton end-point measurement and the model dependence 
may be correlated as well. These correlations will be investigated. This analysis is the 
first to use neutrino reconstruction and the full three-dimensional kinematic information to 
extract the B —>■ Xulv branching fraction. This also pioneers the application of the multidi- 
mensional cuts laid out by Bauer et a/. . A more detailed study of the region of sensitivity 
is in progress, and may result in a less model dependent extraction of \Vub\- 
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The composition of the B — > XJiy is also currently being analyzed and will provide new 
information on the branching fractions, the HQET parameters A and Ai, and the Standard 
Model parameter \Vcb\- 

We gratefully acknowledge the contributions of the CESR staff for providing the lumi- 
nosity and the National Science Foundation and U.S. Department of Energy for supporting 
this research. 
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